function output = mainVlstar(y,x,q,c) 

N=size(x,2);

[betas,gamma,c,fX,e] = estVlstar(y,x,q,c);

output.e=e;
output.betas=betas;

reshapeBetas=reshape(betas,[4 N])';

output.betasRec=reshapeBetas(:,1:2);
output.betasExp=reshapeBetas(:,3:4);

output.gamma=gamma;
output.c=c;
output.fX=fX;

end
